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Ec ABSTRACT 


Parameter plane techniques are applied to compensation of third 
order control systems by extending the principle to include the quad- 
ratic form of the parameters. 

Equations for constant bandwidth curves on the parameter plane 
are developed, 

Recursive equations for the extension of these methods to higher 
ordered systems are presented. 


A computer program for the computation and plotting of constant 


qe 
‘ gt 
6 24 £ 


zeta, omega and bandwidth loci is written, 
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1, Introduction 

Until 1959 when Mitrovic [1] introduced what is now called the 
Coefficient Plane, analysis and compensation of control systems were 
accomplished primarily by frequency response and root locus tech- 
niques. With the advent of the Coefficient Plane there was then a 
method available which allowed the control engineer to determine the 
root locations of a plant as two coefficients of the charactcristic 
equation were varied. 

At about the same time E.R. Ross and T.C. Warren [2] in 
conjunction with Dr. G.J. Thaler published a numerical method of 
deternining single and multiple section cascade compensation while 
maintaining a constant error coefficient, Shortly thereafter in 
1960 C.D. Pollak [3] and Dr. Thaler published yet a different 
numerical method for stating compensator pole locations for a 
given gain, 

In 1964, Siljak [4] extended Mitrovic's work to apply to pro- 
blens in Which the system variables or parameters appeared as a 
product in one or more coefficients. This paper will further extend 
the Parameter Plane concept to include the quadratic form of the 


coefficients. Specifically the case of: 


ч k 
Ж S = 0 (1-1) 
]c=O 
where: 


a, = by A + су. + а, /5 + е,/9 + 1,8" + Gy, (1-2) 
and: 


б, эс,»4, е f 8). аге real and constant 


c and 5 are real and variable 


These curves which will be plotted on the c = 6 or Parameter 
Plane will specify two of the roots of Equation (1-1) for any 


choice of с апа /8 ° 


2, Concepts of the Parameter Plane 

The Parameter Plane was derived by a study of the character- 
istic equations of control systems where two parameters were allowed 
to vary. Here a presentation of the basic ideas of Siljak [4 ] viii 
be studied, 


Consider a characteristic equation of the form: 
n 
2. а в“ = 0 (2-1) 
k=0 
where the coefficients 81» Βργοο»» а, аге real and s is the 


Laplace Transform variable: 


з = С + јо = - Jw + jo) 1- 37 (2-2) 


By forming powers of s it is shown [4 ] that: 


wi rod) + ul} (2-3) 


where: 


»(-3)» (-1)ΐπ( 3) 


(2-4) 
w(-$)- (477565) 
The 9, and T. are the first and second kind of Chebyshev 
functions and are generated by the relations ЕЙ : 
IO 296209) «тү (37”-70 
(2-5) 
with: 
T($) = 1 008) = ο 
p) - 6 (3) -- 


9. 


Taking Equation (2-1), substituting Equation (2-3) in it, then 
equating the real and imaginary parts to zero separately and employ- 
ing Equations (2-4) and (2-5) we will obtain the expressions in 


terms of Chebyshev functions: 


n 


28 (-Иға о (2) = 0 (2-6) 
К-О 
> (а) аа (3) = o (2-1) 


Now, assuming as Siljak did, that the а. coefficients of 
Equation (2-1) are linear functions of A and 9 : 
а = dA + opt a, (2-8) 
with: 
bis O5 6 constant and real 
cx and /В variable and real 


we may then substitute Equation (2-8) into Equations (2-6) and 


(2-7) and obtain: 


BX + pP = -D, 


BLA + с,/8 = -D, (544 
where: 
š y Cop 7 qu acabas (δι. ἐστ 
ο. тов ал κ - k^ κ 
E ης τοι ауу 0, = > ET a 2 
“1 2 À сл 2 2 ία) ο T (2-10) 


а 
li 


II 


1 = (-1)Fa cs U, , D, > (туа ш. 


k=0 


10. 


Now Equations (2-9) may be solved simultaneously for @ апа /3 


giving: 
ons Eu a 
Санса ü / - aa (2-11) 
vo u ұза "%% 


Equations (2-11) now give œ and /8 in terms of 5 and W and 
either may be fixed and the other varied from zero to infinity which 
will give a curve on the οἱ - /8 plane, This may be repeated for 
as many values as needed giving an entire family of curves on the 

сх - 3 plane which completely describe the behavior of a pair 
of complex roots, These curves will usually be constant zeta or 
constant omega loci. 

If s= G is substituted in Equation (2-1) and the same 
procedure is followed, curves of constant sigma or values of s 
along the real axis may be obtained. 

This then is the general procedure to be followed in the devel- 


opment which follows for the quadratic case, 


ll. 


3, Cascade Compensation of Third Order Systems 
The initial investigation of the ox - 8 quadratic case is made 
by attempting to compensate a specific third order system with two 
identical sections of cascade compensators maintaining a constant 
velocity coefficient. The uncompensated plant, Figure 3-1, is chosen 
so as to afford not only lead, but lag compensation, and with this 
selection, the system gain is set to place the complex roots at the 
Stability limit, which also insures the need for double section comp- 
ensation over a large range of zeta and omega, Establishing the 
generalized compensated plant, Figure 3-2, the characteristic 
equation is determined: 
PME τ + 2р2) 84 + (52° + 12p2°) a> + (10pz* + 
“2. + 20 уз + (5p*z* ү бора) в + 30022 = 0 (3-1) 
Simplifying this equation and making a substitution of variables: 
je = p x = p/z (3-2) 
an equation in terms of c and /8 is obtained: 
в? + (6а 2/8 )s* + (5 + 12/3 :/84) а + (10/9 + 63° + 
30«*)s* + (5/37 + 30 о/8 )в + 398° = 0 (3-3) 
and when the terms are arranged in terms of the parameters: 
42 (3085) + о/8 (боз) + pst + 128° + 108°) + 
„В (е η 58 + 30) 4 (s? + 68" 4 582) = 0 (3-4) 
which shows the general coefficient form: 


n 


Ss (ъ + 0 A + 2.2/8 + ex^ + 2/8 E в.) - 0 (3-5) 


k=0 


12. 






K 2 
5(5-1)(5--5) 






UNCOMPENSATED SYSTEM 


FIG. 3-1 


ЫС, P (Sez) τ. 390 
ad 2? (S +P)? 5(5%1)(5%5) 


: COMPENSATED SYSTEM 
FIG. 3-2 


13, 


ог: 
n 
k 
DS sus - ο (3-6) 
с а 
where all of the c's are zero, 


Employing Equations (2-6) and (2-7) 


E» (-1) N a = O (2—6) 


k=0 


>. (-1)%, оз, (5) 


ll 
© 


(2-7) 


Equation (3-4) may be split into a real and imaginary part in terms 


of 1 and ¿Sm 


A Boy + 2/8 Ὁ) + AE, + B Pa “aime ^ . 
AB + 061, + въ + Въ, + = 0 (3-8) 
where: 
By, = 30a°U,( 3) (3-9) 
В,, = 30 a^U,( S ) (3-10) 
D, = -60%0 (3) (3-11) 
D, = -60eU,(3) |. (3212) 
By =20%,($) - 120%,(3) + 100%, ($) (3-13) 
Ej, - 200,3) - 120, (3) + 100 ,($) (3-14) 
1 


The second number on the double subscripted variables and 
the number of the single subscripted variables indicates the real(1) 
or imaginary(2) part of the equation. 


14. 


Fj 
ti 


оү = ($) + 600,(3) - 5wu (3) «3ου (5) (3-15) 


Pop = =w UU ($) + 60 US$) - 500, ($) + 300 (5) (3-16) 


G, = WU (3) + cau (3) - 55) (3-17) 
G, = -wU,($) + owt, ($) - 5г70.(3) (3-18) 
Further collecting the terms of Equations (3-7) and (3-8) we 
obtain: 
A Bos + AP, + Q4 = O (3-19) 
ос “Ву, «ЯР, %9%, = O (3-20) 
where: 
P, = б (8.21): 
ЕЕ - Въ, (3-22) 
ὃ) = E + Чы. + б, (3-23) 
| 22 
q, = ВЕ +f Er Es (3-24) 


We now have obtained a set of two non-linear algebraic equations in 
terms of á and O2 which must be solved simultaneously for the 
correct values of «A ала /9 . То до this, the method with the best 
chance of sucess appears to be Sylvester's Method in which we form 
a set of four equations by taking the original Equations (3-19) and 
(3-20) and forming two more by a multiplication with & giving: 


2 
A BJ] + AP, + Q 


0 (3-19) 


i 
© 


οἱ By, + XP, + Q, (3-20) 


1». 


xB + өсер + AQ, 


21 1 


3 2 
«Во + КР, 


Now placing these equations 


O Boy 


O Boo 


шог Ug 


В. Po 


we may say that if the 


+ X lo 


І 
O 


О 


Pe 
Po Ὁ 
Q. O 
Q5 O 
cx 's are 
T A 
Po ®% 
Q 0 
Qo O 


in matrix form: 


not zero then: 


Expanding this determinant we obtain: 


2 2 
-85145 % 81822814: % Р1Р8,14: - RM + 9192821822 


22 2. 
AS 


2) + «EP 


28552 = 9 


(3-25) 


(2-6) 


(3-27) 


(3-28) 


(3-29) 


Substituting Equations (3-21) through (3-24) in Equation (3-29) we 


achieve a fourth order equation in б : 


16, 


4 2 T тл ту та = ооо 2 - 
A (йз ЖЫЛКЫ жыш Еш л 


LI E 


| 2 


028ү + 01051185) + 


3 n т 2 та p 
/87егв,2 2592 528188182 % 20111281822 + 


2 2 
МВ 5821 7 2: ВъФвоо - “Бъд 2 Boo — 


2 
DSB] Bo] + 2,228,182) 3 


82-87282 _ om Gf l mun n p + 6,8. В 


20722701 Ы 12721722 μμ t 
OG FB B. +DDG3B. = DGB. - E2 p° - 
3500902253 WR 27 " oe 
25 p He ' атг ) + 
оу буро» Ше ән % 21024822 
jJ us 2 
/A(-?25,50,85, * 2E,,058,,855 4 2G,E, 535,8, - 23110182) + 
2 2 2 
(-282) + 20108,18, – 0185) = 0 (3-30) 


from which the coefficients may be determined by a substitution of 





fons (3-9) through (3-19) and the values of the first kind of 
Chebyshev functions, Equation (2-5), in terms of 3 and WO . Since 
the solution of a fourth order equation is at best difficult, it is 
at this point a digital computer becomes a necessity. 

The major problem is not the actual solution of the quartic 
itself, but rather the proper choice of one of the four solutions, 
There are two marked characteristics, however, which help in the 
selection. These are: 

ae Complex answers to the quartic have no physical significance 


and may therefore be discarded as erroneous. 


11. 


b, The definition of GK requires that A and /9 be of 
the same sign so that p and z will be of identical sign. 
Using this information and that available from the Ross—Warren [2] 
method (Appendix I) as to compensator pole and zero location, it is 
found that the solution to the /8 quartic (Appendix II) is the 
largest, positive, real value, 
Now entering Eugation (3-21) with this value, and evaluating the 


other coefficients (Appendix II) we find that: 


οἱ - É Q, / в, | % (3-31) 


for in the third order case Р, is always identically zero, 

Thus, with the programing of the appropriate equations, the 
digital computer could give all of the values and plot the constant 
zeta and constant omega loci on the Parameter Plane for any desired 
values. With the œ and P values, Figure 3-3, and the graphs of 
the Parameter Plane Lag, Figure 3-4, and Lead, l'ígure 3-5, compen- 
sation zones, a spot check of the velocity coefficient is made with: 

2 
22702175 (77: 
V s>+0 с c | P. TT Pi 


where Ko is the compensated root locus gain, 


(3-32) 


Case 1, $ -.5 o» 21.0 Ps —.5 + ј 0.86603 


с 


Д. "НІ? 
К (1.0)(.5* n .Β66ο35}5(4.55 " .866032)2(.469632 + .866032) 
(.41769° + .866032) 


p = ñ - 0.03037 & = p /Zg = 0.36895. 2, = 0.08231 


K, = 4.81100 
2 
бы» (4.81100)(.09231) ) 7.06 
ΠΠ. ΠΡ, 


Ше, 
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Case 2, 3-.2 65-30 г--.6ж 3 2.93943 
/9 = p, = 4.5954 c -p,/2,-1.131 2, = 2.65002 


El. 
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! 
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Ko = 90.21638 
2 
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Case 3, d E .8 ОЎ = 4.5 г 


ИЗ = р, = 10.191 % = p. / 2, 


Зе 


3.373 de = 3.02134 
i. 4. 
К - 5 2.62 = αμ. = 1.4° + 2, й E бъ пе + 2.15 
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с 


K. = 341 „3239 
2 
к, = 341.3239) (3.02134) _ 6,0015 
(2.0) (5.0) (20.299) 


Noting that all of the calculated K's are greater than or equal to 
6.0, the velocity coefficient of the uncompensated system, & check 
of the root loci is made. As shown in Figures 3-6, 3-7 and 3-8, 


the desired root location is also achieved. 
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4. Generalization of the Third Order System 

In Section 3, the development of a specific type one, third 
order system was shown. Here a generalized development of the 
equations necessary to write a computer program for any third 
order SHhten is shown, 

Starting as before by specifying the compensated plant, (this 
time however in generalized terms), as in Figure 4-1, the character- 


istic equation is written. 
в? + (2p. 4E г.) 9 + (ве + гр г. «2 Tr,)o? + 


(в т, + ails, + Tr, )s° + ГрАЩА + гр г.) в + 


2 Е 2 5 2 
pz, + К(Ё5 в + 2Е s + p°) (4-1) 
2 
where: 
EA а the sun of the roots 
Tr, Ф the sum of the product of the roots taken 
Ёл h ata time 
Tr. 2 the product of the roots 


4. 


With the substitution of: 


f= p © = р/2 (4-2) 


and the collection of like terms in eX and /8 , не obtain: 
2 2 2, 3 2 
ο (Кз) + a /S(2Ks) + © (s + г.в +2, Trr,s + Пт. + К) 
4 3 2 
+ Αί25 + 20.2, в + 22 Tr,s + 2Tr, s) 


+ (s? + Bris tA Trs a Trs’ (4-3) 
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Employing Equations (2-6) and (2—7) and defining: 


В, = Кош (3) (4-4) 
в,, = Kuu ($) (4-5) 
Dj, = -2Kau,($) (4-6) 
D, - -2Кой (3) (4-1) 


E] = 2 «0. (S) -22 7,0%, ($) + 2L.Tr, wt, ($) 


- 2ilr,0u,($ ) (4-8) 
Ер“ 2 ο (5) -20,г,0970.(3) + 22, s uo? ($ ) 

- 212, 00.03) (4-9) 
ғ --о%3)- Аты(3)-2Тт,өз, 2) 

+ (Те, + ко. (3) (4-10) 
Po = US) 27,03) -Z Mron ($) 

+ (Mr, + K)u_($) (4-11) 
G = -w ($) +Z z; oot ($) – ТГ, о 20,03) 

+ Tr wu ($) (4-12) 
в CH rt IT. 

+ Tr w? S) (4-15) 
Ру = Въ (4-14) 
P, = AD, (4-15) 
Q = B®, 4 Bo + б, (4-16) 
Q. = /882% АВ тор ы: (4-17) 
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we find: 
2 
«Bj, * XP, *4Qj = 0 (4-18) 


2 
AB, + AP, + Q, = O (4-19) 


which are two non-linear algebraic equations completely general- 
ized in terms of the uncompensated system poles and root locus gain, 

5 ‚ & and the first kind of Chebyshev Functions. Equations 
(3-30) and (3-31) are now employed to obtain the solutions for A 
and /3 ο using the same selection rules as before. 

Several third order systems were investigated by the application 
of the generalized equations and the Parameter Plane curves, Figures 
4-2 through 4-9, were plotted. ОҒ these, the к/з? family appears 
the most interesting. Further investigation of three of the curves 
in this family, Figures 4-2, 4-3 and 4-4 shows that there is a rela- 
tionship between K, the root locus gain, Οἱ and „В : 

These relations are: 

a. Choose a point on the 1/83 о -/3 plane, 

b, Zeta reads directly. 

c, Determine the actual omega at that point by multiplying 

the value read by the cube root of the uncompensated 
system gain. 

d..Read the value of A directly from the point chosen. 

e, Read the value of 2 from the point chosen, 

f, Obtain the true value of 8 by multiplying this value by 

the cube root of the uncompensated system gain. 


By this method, the values of & and „ИЗ пау be determined for 
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systems from the one universal curve, 

Although the majority of the curyes are easily interpreted, 
there are some, Figures 4-9 and 4-10, which require careful corre- 
lation of the data points and the curves. As shown in Figures 4-11 


and 4-12, the correct curves may be picked out and redrawn if 


necessary for further clarity and ease of use. 
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5. Bandwidth Curves on the аә - / Plane 

In many instances, there is also a bandwidth criterion imposed 
on the engineer as well as an optimal operating point for the plant 
under consideration, With this in mind, equations for the plotting 
of constant bencwidth curves on the oœ = 78 plane are developed. 
For the purpose of this development a constant bandwidth curve will 
be defined as: 


A constant bandwidth curve for G(jw,) - lis & curve drawn 
upon the parameter plane which specifies the relation between 
the paremetors necessary if the transfer function G(s), which 
is a function of the parameters, is to have magnitude M at 
the real frequency 0) р. [ 
Once these curves are obtained they may be superimposed on the 
parameter plane thus indicating what values of the parameters are 


necessary in order to meet the specifications. 


Taking the rational transfer function and defining it: 


ps 4p ae ον p45 * p 

P(s m- пай un ο 

G(s) = пра . E E ЫС (5-1) 
4,5 + а „13 + 000958 + q. 


where the Ра 8 and 4,18 are of the form: 


2 : 2 
р, 6,9 + h o + ip + 21.06 + 668 + 1, 


М = 0,1,2, „еў (5-2) 
2 2 
q. = a. + Dx + caf + 8 В + ey ++ 
v = ο. ὃς ος 9 Ὦ (5-3) 
(е. m d 
we may write: ра pS 


G2 
ms 
2 
< 

l 


> " (5-4) 
W me 
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Employing Баџафіоп (5-4) we write in the parameterized form the 


generalized compensated third order transfer function: 


acie et (5-5) 


с 252 + 29/3 Кз + BK (5-6) 


where: 


P(s) 
ands 
2e 2 
Q(s) - /9 [s +2 Wr,s “С.Т т,» * 220 0 г. | + 
Аре + 22 To, οὖ + 22 Tr, s + LT r,s] + 
52 4 & Mr, VT s.s] + 2 Пее (5-1) 


Making the definitions: 
= EV, V 
2 (-1) со a 5 eto. for B,,C .,D ,E .,F. (5-8) 


even 


m 
lI 


= %(у-1),,У 
A. = Σ, i-i Was ete, for B,,0,,D, ,E,,F; (5-9) 
odd 
= žu u 
С = Σ, (-1)""w g, š eto, for H ,I J. K ,L. (5-10) 
u= 
even 
> 3(u-1) „u 
G. = Zu (-1) 3, 8,3 etc. for H,,I,,J,,K,,L, (5-11) 
odd 


and substituting in Equation (5-5), we write: 
2 
(x G, + K.) + )е/8т.) 
(820 + BE + г) + 8“. + дк. + е.) 
ο 208 2; ΑΝ, i 


Setting the magnitude of с(јоо,) = М ме Ріпа: 





G( jw, ) Е (5-12) 


43, 


2 л? 
а A аа 
(4%, + ВЕ. + г.) + 7:52 + BE, + г.) 


Manipulating Equetion (5-13) algebraically we find that: 


ΝΕ = КЕСЕ = 


ф (4,8) - i (5,8) = 0 (5-14) 
where: 
pap) = x “G° 4 2x K G, + К + 058 212 SN 
9 (οἱ) - B D + 2/87D,B, + 2 8° n, «Απ» 
283 F, + p. + r + BEL + + 
2 8%,D, К Во; + 2 BEF, (5-16) 


Substituting Equations (5-15) and (5-16) in Equation (5-14) and 


defining: 
Р, = ре х Dé (5-17) 
Q, = 20 Е + 2E,D, (5-18) 
R. = 2D P + EZ + E; + 2D;P, (5-19) 
ТЭЭР T L TL 
Ше «ес (5-20) 
V, в 26 Е + 2E,T, (5-21) 
W = Ес + F (5-22) 
Wy =s «Λο καν (5-23) 


we can write: 


44, 


Къв " iae + (ИВ. — RB.) fe + кту8 + 
(мн, — W.) = O (5-24) 

Since the Parameter Plane for compensation purposes has already 
been determined it is now a matter of taking the computed & values 
and substituting them along with a constant value of omega and M 
into Equation (5-24) and then solving the (8 quartic. This has as 
its solution the largest, real and positive value of the four roots 
as before. 

Examples of these curves, drawn so as to enable their use as 
overlays on Figures 4-2 through 4-10 are shown in Figures 5-1 


through 5-9. 
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6, Programing of the Compensation and Bandwidth Curves 

Appendices III and IV are listings of the programs developed 
in the course of investigation of this work. Only PROGRAM PROJECT 
(Appendix III) will be explained, as PROGRAM PLOT (Appendix IV) is 
essentially PROGRAM ROOT LOCUS [6] with a minor subhead and internal 
modification. 

PROGRAM PROJECT is designed to solve the οί quadratic and 
quartic established in Sections 4 and 5. ‘The program is divided 
into two main sections, the first for the computation of the οἱ -6 
points and the second for the bandwidth points, 

The first section computes an 80 by 80 matrix of the A and 
8 points corresponding to set values of Š and € „ The computa- 
tional part is followed by two distinct graphing sections, one for 
lag and the other for lead compensation, 

The lag graphing section is set up so that during the plotting 
of the curves each value of A is tested to determine if its value 
is 107 5 ЄС 1001 , If mo points are found within this range 
then a print out is made: 

NO LAG COMPENSATION POSSIBLE 

For the lead section graphs, @ is again tested by the criterion 
1.00015 A £ (X-graph scale)(X graph width) . Again if there are no 
values of CX within this region the statement: 

NO LEAD COMPENSATION POSSIBLE 
is printed. Іп this case however a study of the printed values of 
Сі must be made to insure that the points are indeed non-existant 
or rather just lie outside the range of the graph. 


The second main section of the program computes the value of 
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/8 for a given value of A as determined by the X graph scale, 
Here the plotting routine is set up so as to not plot zero points 
and to stop the curve when either the A or /8 value exceeds the 
range of the graph. 

There are several go-no go options available within the 
program which are amply explained in the comment section preceeding 


the program listing in Appendix ІІІ, 
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T. Extensions to Higher Order Systems 

Although the work presented to this point has been limited to 
third order systems and the program written for this specific case, 
investigation shows that generalized equations may be written which 
will allow the extension of the program to higher ordered systems. 
It can be shown that for a given ni order system with no zeros to 
be compensated with two identical sections of cascade compensation, 
that the characteristic equation of the system may bc generally 


written as: 


ВМ аш £ ν᾽ 144. 2pzs + Se n 
j=n 1-0 
2l il 2 К-О k 
2 21 (21: )8 «р 21 (Ar) + 
Кеп 9 =n—1 2 
j= 1-1 
j=n 
к-г 
> (Т2, ) з" = Ὁ (1-1) 
К=п+1 д 
jel 


where for n=4 the equation would be written: 


3° 4 ор? + prat casu 2pzs 4 p^)K + 


ar(LT r,s" a 2 2,5 “АТ r,s“ +2 17 г.в) + 
MN ERES к Тв РАН, 272) + 
κ... АТ 7, s^ + Те” XT 2s) - d (1-2) 


It may be further shown that the parameters defined by Equations 


(4-4) through (4-13) may be written: 


2 
Bo] = Ко (3) (155) 
Boo = козу, ( 3) (7-4) 
D, = -2K0U (3) (1-5) 
D, = -2KwU, ($) (7-6) 
ыг 
ο ο ση 
j=1 (7-7) 
jen 
m+] ntl $ p k ky ў 
взе = 26-0) Een ERNEST Те)! 
3-1 (7-8) 


j=n 
k=0 

ш = (-) ауу, 109 ) ЗА. [CASA $ ΣΤ.) ко (S ) 
3-1 (7-9) 


jen 
k=0 

E ID E VADER] as 
k=n=1 d | 
j=1 (1-10) 


j=n 
k=2 
a AT ($)« 2 | (17565, „(3 ET) 
3-1 (1-11) 
j=n 
k=2 
+ Sis : 
De 2 [ao (CTs) 


j= (7-12) 


n+2  n+ 
ы-ы шз", а 
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These then are the recursive equations required for the complete 
generalization to a lan order system. By employing the above 
equations and replacing in PROGRAH PROJECT (Appendix III) cards 100 
through 150 and 300 through 540 with the appropriate programing, the 
program may be used for any given nod order system. 


In like manner by generally defining: 


P(s) e» οι κε; + ос / Кв + 82 K (7218) 
ала: Е: 
Q(s) = uc 2B + 8 + 12 Ë: ЛЕЗДЕ 
Je 
can j=n 
k=2 K 
в Е QW Із; pee D | (214) (1-14) 
1-1 3-1 


and using Equations (5-8) through (5-11) we may replace in the 
program cards 2860 and 2880 through 2920, thus adapting this part 


of the program to a general Қы отдет application, 


8, Conclusions 
“arougnout tnis development of the Parameter Plane quadratic 
extension, the οι.’ 8 in the generalized coefficient form: 


$^ (ъс + са + de 4 efe + г В +g) = O (8-1) 


κο 


have been identically zero. This at first appearance might seem to 
detract from the generalization. The inclusion of this parameter 
does not however introduce any great difficulty in the solution. 
The change in the development would be to the value of P; and Р, 


which would become: 


P. = C1 E Въ (6-2) 
Р, = Cy + Въ (без) 
and the final solution for «A which would change to: 
Р Р 22820210 
- 28 
.--2- 00105. (8-4) 
2B 2 
21 4821 


For this case, new selection rules for acceptable values of cx 

would be used, and would be much like those presented for (8 ° 
Though the extension of the Parameter Plane to include the 

α - /8 quadratic case makes this tool even more useful, further 

work is still to be done in this field. Not only must the equations 


for the solutions of the Parameter Plane curves for such cases as: 


a = b, e *9* + оро 76 + mt" + e, * + t8 + 
6.8 + x + uS юм [5] (8-5) 
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and higher ordered combinations of the parameters be developed, but 
more efficient programing techniques must be developed, In the use 

of PROGRAM PROJECT, for instance, as the location of the system poles 
on the QO axis of the S-plane move to the left, the computational 
time becomes excessive due to present programing technique and computer 
speed. 

Another major problem in further extensions of these techniques, 
and indeed even other applications of the curves from the preceeding 
development, will be interpretation, Іп this case, the initial sub- 
stitution of variables immediately allowed interpretation of the 
curves sight unseen, Here then, will be most likely the one single 
drawback to further extension, for as the parameters & and /9 аге 
used as representations of other variables in control systems, each 


application will have its own unique interpretation. 


6l, 


2. 


3. 


4, 
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APPENDIX I 
THE ROSS-WARREN METHOD OF COMPENSATION 
The Ross-Warren [2] Method is an algebraic way of determining 
the location of compensator poles and zeros for single or multiple 
sections of cascade lead or lag compensation, The 1: нь 
ОҒ the method are based on relationships between the root locus gain 


and the error coefficients, where the uncompensated root locus gain 





n 
: n 
S | TP |z + Р; | 
Ix] - lim (1-1) 
ТГ |r + 2; | 
i=l 
and the compensated root locus gain at the point r is: 
k 
K (2) = K(r) (I-2) 


where К is the number of identical sections introduced by the 
compensator. 
The other part of the equality for the derivation is determined 


by the fact that the uncompensated error coefficient is defined by: 


Js 





i=] а 
K (z) = Қт) (1-3) 
ATUS Mer 
il 1 
and the compensated error coefficient is: 
p Z 
d с 
ШЕ) = Kan) z p. (1-4) 


These relationships are expressed geometrically on the Splane in 


Figure I-l and are used to derive the relationships: 
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S- PLANE RELATIONSHIPS ТО 
DERIVE THE THE ROSS- WARREN EQUATIONS 


<- 0, + jw, 
κ. 
3 x 
1 
ж 
292 о 
\ 
\ Ө 
A : 3 N „е 
р 2, 
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а = al sin A (1-5) 
67 Қ(т) sin © = 


h sin À 


Ро 7 sin (8-4) (1-6) 


where: 

λος; (οοἱ φ- τπτ ese ᾧ ) (1-7) 
and: 

Ө-Т»%-3-0 (1-8) 
Неге: 


h ë = the distance from the origin to the root point 
and is specified by the choice of zeta and omega 


K = root locus gain of the uncompensated system 


K(r) - root locus gain of the uncompensated system at 
the point r 


ф - angle required to make r a point on the root locus 


A - an angle calculable from K, K(r) and ф 
3 - angle between the negative real axis and h 
Ө – an angle calculable from Q , 3 and A 


These equations are generalized for n sections so that: 


d sin À 
K i^ d 
| у sme син 
h sin A 
n 


where: n 
An s qa" (cot Φ., - үсу ЕЕ ф P (I-11) 


теъ" ф,- б = у (1-12) 
апа: ф / a m 


Ф 
O 5 
2. II 
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APPENDIX II 
DETERMINATION OF BETA QUARTIC SOLUTION 
Here, three examples of the approximately fifty calculations 
that were made in the investigation of the proper solutions of Οἱ 
ала /9 » will be presented, "hose calculations are started by em- 


ploying Equations (I-9) through (I-13) with n equal to two. 


Case 1, ENS 10843 1.28 
$ - tan | =. ten | 2:5 + tan” #5 ч 1809 = 24938! 
O 
Ф, = 12919: f 
K(r) = (4.0)(17.0)2(41,0)2 = 105.60304 
[κ/κ(π) = .532994 
A> = cot" (cot 12219! - (.532994)ese 129191) = 25040! 
Ə, = 180° + 120191 - 90% 25 доп та 462301 
4 ο 
p. =Q sin 250. gaps 
sin 64°20! 
Z = EG 532994) in 25 40! - .94906 
sim 76° 30! 


this implies that: 


6 - 1.9225 о6 - 2.02568 
Case 2, oF = Gar Ke 30.0 h = 5.0 
Ф = tan (453202) y dam (423302) y tam (453302) _ 180° 
b = 10916! 
O 
$, = 54 33! 


- нЕ: 
К(г) - (5.0)(21.000199)%(25.00019)% - 114.8164 


lx/x(x)" = «511029 
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ЩЕ = соі! (со% 54733! - (.511029)сзс 54933') = 85°09.5! 
E ТО 5254? 33: — 7699 — 85809,5: - "ВО 2.5: 
e 1 > е 
p. = Gre. sin 85 09.5! - 8.12051 
sin 34 ο 
= (5.0)(.511029) Mi μα. - 2.54616 
sin 89? 23,5! 
this implies that: 
| /5 = 8.72057 A = 3.42498 
Case 3. 31-00 Bae 3070 ПЕЕ 5.8 


tan" Қа + tan 1 (2) + tan es Q^" e ppt 


$ 


ф, - 1654 
d д. 
K(r) (560)02640) "(5080)" = 1808277586 


YK/K(r) = .407934 


N, = cot” (cot 16%541 - (.407934)ose 169541) = 27°59! 
Ө, = 180" «16754! - 90" - 2179591 = 1855 


ο 
р. = „О 225 21 ! 2 2.60005 
sin 62 01! 
b —sin 21594 
2. = (5.0) („дотоза)--В Е = 95444 
sin 78° 55! 


this implies thats: 
/Ӧ = 2.60005 & = 2.72416 
Now taking the values of /8 just found we compare them to the 
roots of the /8 quartic given in Figures II-lb, II-2b and II-3b 
respectively. We see that the proper answer for /? y in all cases, 


is indeed the largest positive and real root of the quartic. From 


Si. 


Figures II-4 and II-5 we now obtain the values necessary to solve 
for the value of A . These steps for the three cases are shown in 
Figure II-6, and again we note that only when applying the above 


selected values of (à do we obtain the correct value of CX , 
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